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Propagation of light in the gravitational field of self-gravitating spinning bodies moving with arbi- 
trary velocities is discussed. The gravitational field is assumed to be "weak" everywhere. Equations 
of motion of a light ray are solved in the first post-Minkowskian approximation that is linear with 
respect to the universal gravitational constant G. We do not restrict ourselves with the approxima- 
tion of gravitational lens so that the solution of light geodesies is applicable for arbitrary locations 
of source of light and observer. This formalism is applied for studying corrections to the Shapiro 
time delay in binary pulsars caused by the rotation of pulsar and its companion. We also derive the 
correction to the light deflection angle caused by rotation of gravitating bodies in the solar system 
(Sun, planets) or a gravitational lens. The gravitational shift of frequency due to the combined 
translational and rotational motions of light-ray-deflecting bodies is analyzed as well. We give a 
general derivation of the formula describing the relativistic rotation of the plane of polarization of 
electromagnetic waves (Skrotskii effect). This formula is valid for arbitrary translational and rota- 
tional motion of gravitating bodies and greatly extends the results of previous researchers. Finally, 
we discuss the Skrotskii effect for gravitational waves emitted by localized sources such as a binary 
system. The theoretical results of this paper can be applied for studying various relativistic effects 
in microarcsecond space astrometry and developing corresponding algorithms for data processing in 
space astrometric missions such as FAME, SIM, and GAIA. 



I. INTRODUCTION 



The influence of gravitation on the propagation of electromagnetic rays has been treated by many authors since 
Einstein first calculated the relativistic deflection of light by a spherically symmetric mass Currently there is 
much interest in several space missions dedicated to measuring astrometric positions, parallaxes, and proper motions 
of stars and quasars with an accuracy approaching 1 microarcsecond (/xarcsec) Q . Thus progress in observational 
techniques has made it necessary to take into account the fact that electromagnetic rays are deflected not only 
by the monopole gravitoelectric field of light-deflecting bodies but also their spin-gravitomagnetic and quadrupolar 
gravitoelectric fields. Propagation of light in the stationary gravitational fields of rotating oblate bodies with their 
centers of mass at rest is well known (see, e.g., || - Q, and references therein). Quite recently significant progress 
has been achieved in solving the problem of propagation of light rays in the field of an isolated source (like a binary 
star) that has a time-dependent quadrupole moment and emits quadrupolar gravitational waves j?]]. However, the 
translational motion of the source of gravitational waves was not taken into account in Q); in effect, the source was 
assumed to be at rest. 

We would like to emphasize the point that most gravitational sources are not in general static and move with respect 
to the observer in a variety of ways. As revealed by approximate estimates the precision of planned astrometric 
space missions 0, binary pulsar timing tests of general relativity very long baseline interferometry |lC| ] - [CD, 
etc., necessitates the development of more general methods of integration of equations of light propagation that could 
account for such translational motions. Indeed, non-stationary sources emit gravitational waves that weakly perturb 
the propagation of electromagnetic signals. The cumulative effect of such gravitational waves may be quite large 
and, in principle, could be detectable fl2|| . It is necessary to know how large this influence is and whether one can 
neglect it or not in current and/or planned astrometric observations and experimental tests of general relativity. 
When considering propagation of light in the field of moving bodies it is also worth keeping in mind that gravitational 
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interaction propagates with the speed of light in linearized general relativity [[l3| - ]Ti}| , This retarded interaction 
has important consequences and can play a crucial role in the theoretical prediction of secondary relativistic effects 
in time delay, light deflection, polarization of light, etc. 

The first crucial step in solving the problem of propagation of electromagnetic rays in the retarded gravitational 
field of arbitrary- moving bodies has been taken recently in ]15[ ] . The new formalism allows one to obtain a detailed 
description of the light ray trajectory for unrestricted locations of the source of light and observer and to make 
unambiguous predictions for possible relativistic effects. An important feature of the formalism is that it is based on 



the post-Minkowskian solution [|16| - 19 of the linearized Einstein field equations. Thus the amplitude of gravitational 
potentials is assumed to be small compared to unity, but there are no a priori restrictions on velocities, accelerations, 
etc., of the light-deflecting bodies. In this way the retarded character of the gravitational field is taken into account 



in the linear approximation. This is in contrast to the post-Newtonian approximation scheme |20j - |23|, which 
assumes that velocities of light-deflecting bodies must be small with respect to the speed of light. Such a treatment 
destroys the causal character of a gravitational null cone and makes the gravitational interaction appear to propagate 
instantaneously at each step of the iteration procedure. As proved in Jl5| , this is the reason why the post-Newtonian 
metric gives correct answers for the time delay and light deflection angle only for a very restricted number of physical 
situations. As a rule, more subtle (secondary) effects in the propagation of light rays in time-dependent gravitational 
fields are not genuinely covered in the post-Newtonian approach, whereas the post-Minkowskian light-propagation 
formalism gives unique and unambiguous answers. 

In the present paper, we extend the light-propagation formalism developed in p"5| ] in order to be able to include the 
relativistic effects related to the gravitomagnctic field produced by the translational velocity-dependent terms in the 
metric tensor as well as spin-dependent terms due to light-deflecting bodies. We shall start from the consideration of 
the energy- momentum tensor of moving and spinning bodies (section 2). Then, we solve the Einstein field equations 
in terms of the retarded Lienard-Wiechert potentials (section 3), describe the light-ray trajectory in the field of these 
potentials (section 4), and then calculate the deflection angle and time delay in the propagation of electromagnetic 
rays through the system of arbitrary- moving and spinning particles (section 5). We pay particular attention to the 
calculation of gravitomagnetic effects in ray propagation in pulsar timing and astrometry (section 6). Moreover, we 
discuss the relativitic effect of the rotation of the plane of polarization of electromagnetic rays (Skrotskii effect |24|] 

pq ) in the gravitomagnetic field of the above-mentioned system of massive bodies (section 7). This effect may be 
important for a proper interpretation of a number of astrophysical phenomena that take place in the accretion process 
of X-ray binaries (2^] - |2^] and/or supermassive black holes that may exist in active galactic nuclei. Finally, we derive 
an exact expression for the rotation of the plane of polarization of light caused by the quadrupolar gravitational waves 
emitted by localized sources (section 8). The treatment of the Skrotskii effect given in the final section may be 
important for understanding the effects of cosmological gravitational waves on the anisotropy of cosmic microwave 
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background (CMB) radiation at different scales |30| - Details related to the calculation of integrals along the 
light-ray trajectory are relegated to the appendices. 



II. ENERGY-MOMENTUM TENSOR OF SPINNING BODY 

We consider an ensemble of N self-gravitating bodies possessing mass and spin; higher-order multipole moments 
are neglected for the sake of simplicity; for their influence on light propagation via the mathematical techniques of 
the present paper see [p3| . Propagation of light in the gravitational field of arbitrary-moving point-like masses has 
been studied in |]l5| . The energy- momentum tensor T a @ of a spinning body is given 

T^(t,x)=r^(t,x)+Tf (t,x), (1) 

where T^f and Tg ^ are pieces of the tensor generated by respectively the mass and spin of the body, and t and x are 
coordinate time and spatial coordinates of the underlying inertial coordinate system. In the case of several spinning 
bodies the total tensor of energy-momentum is a linear sum of tensors of the form (|l|) corresponding to each body. 
Therefore, in the linear approximation under consideration in this paper, the net gravitational field is simply a linear 
superposition of the field due to individual bodies. 

In equation (Q) , T^f and Tg 13 are defined in terms of the Dirac 5- function |34j - |3(| as follows |}7j - ]39[ | 

/+oo 
p^\-g)-^8{t - z°(rj)) tf(x - m(rf))dfi , (2) 
-OO 

r+oo 

Tf(t,x)=- V J S^u^i-gr^Sit-z^v^S^-z^dn, (3) 



— OO 



where r\ is the proper time along the world-line of the body's center of mass, z(rj) are spatial coordinates of the body's 
center of mass at proper time 77, u a (rj) — v l ) is the 4-velocity of the body, u°(rf) = (1 — v 2 ) -1 / 2 , (v % ) = v(rf) is 
the 3-velocity of the body in space, p a {ii) is the body's linear momentum (in the approximation neglecting rotation 
of bodies p a {rf) — mu a , where m is the invariant mass of the body), S a ^{ri) is an antisymmetric tensor representing 
the body's spin angular momentum attached to the body's center of mass, V7 denotes covariant differentiation with 
respect to metric tensor g a p, and g = det(g Q ^) is the determinant of the metric tensor. 

The definition of S a ^ is arbitrary up to the choice of a spin supplementary condition that is chosen as follows 

S a % = ; (4) 

this constraint is consistent with neglecting the internal structure of the bodies involved, so that we deal in effect 
with point-like spinning particles @ - Q . We introduce a spin vector S a which is related to the spin tensor by the 
relation 

S af3 = r] a ^ s u 7 Ss , (5) 



where f] a ^ S is the Levi-Civita tensor related to the completely antisymmetric Minkowskian tensor e a f3 7 s 0] as follows 

V aPlS = -(-gy 1/2 e a( 3 7 s , V a p-yS = (-9) 1/2 e a07 s , (6) 

where £0123 = +1- The spin vector is orthogonal to u a by definition so that the identity 

S a u a = (7) 



is always valid; this fixes the one remaining degree of freedom S°. Hence, the 4- vector S a has only three independent 
spatial components. 

The dynamical law y^T"' 3 = 0, applied to equations (|l|) -@, leads to the equations of motion of a spinning 
particle in a gravitational field. Indeed, using a theorem of L. Schwartz that a distribution with a simple point as 
support is a linear combination of Dirac's delta function and a finite number of its derivatives, one can develop the 
theory of the motion of point-like test bodies with multipole moments in general relativity E3| . In this way, it can be 
shown in particular that for a "pole-dipole" particle under consideration in this paper, one is led uniquely p3] to the 
Mathisson-Papapetrou equations with the Pirani supplementary condition (^) . 

In what follows, we focus only on effects that are produced by the spin and are linear with respect to the spin and 
Newton's gravitational constant (the first post-Minkowskian approximation) in the underlying asymptotically inertial 
global coordinate system. The effects produced by the usual point-particle piece of the energy-momentum tensor have 
already been studied in Jl5[ . Let us in what follows denote the components of the spin vector in the frame comoving 
with the body as J a = (0, J7"). In this frame the temporal component of the spin vector vanishes as a consequence of 
and after making a Lorentz transformation from the comoving frame to the underlying "inertial" frame, we have 
in the post-Minkowskian approximation 

S°= 7 (v.7), S l = J l + ^—}-{vJ)v\ (8) 

where 7 = (1 — t> 2 ) -1 / 2 and v = (V) is the velocity of the body with respect to the frame at rest. 

III. GRAVITATIONAL FIELD EQUATIONS AND METRIC TENSOR 

The metric tensor in the linear approximation can be written as 

g a p(t,x) =T] a i3 +h a p(t,x) , (9) 

where r/afj — diag(— 1, +1, +1, +1) is the Minkowski metric of flat space-time and the metric perturbation h a /3(t,x) 
is a function of time and spatial coordinates p(| . We split the metric perturbation into two pieces h 1 ^ and h'g 3 that 
are linearly independent in the first post-Minkowskian approximation; that is, 

h^ = h^ + hf. (10) 
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Thus, the solution for the each piece can be found from the Einstein field equations with the corresponding energy- 
momentum tensor. 

The point-particle piece h^S of the metric tensor has already been discussed in p5[ and is given by 



(ii) 



r(s) — v(s) • r(s) ' 

where r(s) = |r(s)|, r(s) = x — z(s), and both coordinates z and velocity v of the body are assumed to be time 
dependent and calculated at the retarded moment of time s defined by the light-cone equation 



s + |x — z(s)\ = t , 



(12) 



that has a vertex at the space-time point (t, x) and describes the propagation of gravitational field H] on the 
unperturbed Minkowski space-time. The solution of this equation gives the retarded time s as a function of coordinate 
time t and spatial coordinates x, that is s = s(t,x). 

As for hg , it can be found by solving the field equations that are given in the first post-Minkowskian approximation 
and in the harmonic gauge |l5| as follows 

Uhf{t,x) = -16rr 

Taking into account the orthogonality condition (H), the field equations (tL3|) assume the form 



I^(t,x)--^T| A (*,x) 



(13) 



ST^irfaflfjiMt - z°(r))) <J(x - zfa)) 



(14) 



where we have replaced the covariant derivative V 7 with a simple partial derivative d 1 — d/dx 1 and u a — 7(1, u l ), 
where 7 = (1 — u 2 ) -1 / 2 is the Lorentz factor. The solution of these equations is given by the Lienard-Wiechert 
spin-dependent potentials 



hf(t,x) 



(15) 



_r(s) — v(s) • r(s) 

If we restricted ourselves with the linear approximation of general relativity, the sources under consideration here 
would have to be treated as a collection of free non-interacting spinning test particles each moving with arbitrary 
constant speed with its spin axis pointing in an arbitrary fixed direction. This simply follows from the Mathisson- 
Papapetrou equations in the underlying inertial coordinate system. Thus, in such a case carrying out the differentiation 



in (15), we arrive at 



fcf(i,x)=4(l-*; 2 ) 



(16) 



[r(s) — v • r(s)] 3 ' 

where v and S a @ are treated as constants and we define r a = (r, r). However, we may consider (0) as expressing the 
first two terms in a post-Minkowskian perturbation series; that is, at some "initial" time we turn on the gravitational 
interaction between the particles and keep track of terms in powers of the gravitational constant only, without making 



a Taylor expansion with respect to the ratio of the magnitudes of the characteristic velocities of bodies to the speed 
of light. The development of such a perturbative scheme involves many specific difficulties as discussed in [|l6| - [ ^9[ . 
In this approach one may relax the restrictions on the body's velocity and spin and think of v and S a @ in ( p~6| ) as 
arbitrary functions of time. In this paper as well as p"5| , we limit our considerations to the first-order equations (|lC|), 
(pd|), and (|l6|); however, in the application of these results to the problem of ray propagation, we let v and S al3 in 
the final results based on ( |TT[ ) and ( |l6| ) be time dependent as required in the specific astrophysical situation under 
consideration. Therefore, the consistency of the final physical results with the requirements of our approximation 
scheme must be checked in every instance. 



The metric tensor given by equations (|10[), (|ll|), and (16) can be used to solve the problem of propagation of light 



rays in the gravitational field of arbitrary- moving and spinning bodies. 

IV. PROPAGATION LAWS FOR ELECTROMAGNETIC RADIATION 

The general formalism describing the behavior of electromagnetic radiation in an arbitrary gravitational field is well 
known p6| . A high-frequency electromagnetic wave is defined as an approximate solution of the Maxwell equations 
of the form 

F a p = Re{A a p exp(i(^)} , (17) 

where A a $ is a slowly varying function of position and (p is a rapidly varying phase of the electromagnetic wave. From 
equation (jllj) and Maxwell's equations one derives the following results (for more details see [[l7| - and Appendix 
A). The electromagnetic wave vector l a = dip/dx a is real and null, that is g a ^l a l/3 — 0. The curves x a = x a (X) that 
have l a = dx a /d\ as a tangent vector are null geodesies orthogonal to the surfaces of constant electromagnetic phase 
ip. The null vector l a is parallel transported along itself according to the null geodesic equation 

l p Vpl a = . (18) 



+ r^/ J r = o, (19) 



The equation of the parallel transport ( |l8| ) can be expressed as 

dl' 
~d\ 

where A is an affine parameter. The electromagnetic field tensor F a p is a null field satisfying F a pl^ = whose 
propagation law in an arbitrary empty space-time is 

D dx a 

D\F a fi + 9F a p = , Dx =D\ = ~d\ Va ' (20) 

where 9 = iV Q P is the expansion of the null congruence l a . 

Let us now construct a null tetrad (l a ,n a ,m a ,rn a ), where the bar indicates complex conjugation, n a l a = —1 
and m a rn a = +1 are the only nonvanishing products among the tetrad vectors. Then the electromagnetic tensor 
F a p = Ref^c^) can be written (see, e.g., - ||(|) as 
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= $ l[a m P] + * ka m P] > ( 21 ) 

where $ and \P arc complex scalar functions. In the rest frame of an observer with 4- velocity u a the components of 
the electric and magnetic field vectors are defined as E a = F a ^u^ and H a = (— l/2)e a ^ s F 1 ^up 1 respectively. 

In what follows it is useful to introduce a local orthonormal reference frame based on a restricted set of observers 
that all see the electromagnetic wave traveling in the +z direction; i.e., the observers use a tetrad frame e"^ such 
that 

e Q (0) = u a , e Q (3) = Ha^r 1 [1° + (hu")u a ] , (22) 

and e a (iy e%) are two unit spacelike vectors orthogonal to each other as well as to both e°L-< and e"^ |5lj] - p2[ . 
The vectors e"^ and e" 2 j play a significant role in the discussion of polarized radiation. In fact, the connection 
between the null tetrad and the frame e"^ is given by l a = — (7 7 u 7 )(e" ^ + e a ^) , n a — — i(/ 7 w 7 )(e" Q ^ — e" 3 )) , 
m a = 2~ 1 / 2 (e" 1 ) + anc ^ ™° = \ e °[x) ~ * e "2)) • The vectors e"^ and e" 2 ) are defined up to an arbitrary 

rotation in space; their more specific definitions will be given in section 7. 

Vectors of the null tetrad (l a ,n a ,m a ,rn a ) and those of e"^ are parallel transported along the null geodesies. 
Thus, from the definition (|2l]) and equation (^) it follows that the amplitude of the electromagnetic wave propagates 
according to the law 

— +6»$ = 0, — + 6>* = . 23 

dX dX 

If A is the area of the cross section of a congruence of light rays, then 

^ = 29A. (24) 

Thus, A \Q\ 2 and A \^\ 2 remain constant along the congruence of light rays l a . 

The null tetrad frame (l a , n a , m a ,m a ) and the associated orthonormal tetrad frame e a ,ps are not unique. In the 
JWKB (or geometric optics) approximation, the null ray follows a geodesic with tangent vector l a = dx a /dX, where A 
is an affine parameter. This parameter is defined up to a linear transformation A' = A/A+const., where A is a nonzero 
constant, so that l' a = dx a /dX' = A l a . The null tetrad is then (l' a , n' a , m Q , m a ) with n' a = A _1 n a ; this affine 
transformation leaves the associated tetrad e"^ unchanged. Moreover, we note that at each event, the tetrad frame 
e"^ is defined up to a Lorentz transformation. We are interested, however, in a subgroup of the Lorentz group that 
leaves l a invariant, i.e. A"jZ^ = l a . This subgroup, which is the little group of the null vector l a , is isomorphic to the 
Euclidean group in the plane. This consists of translations plus a rotation. The translation part is a two-dimentional 
Abelian subgroup given by 

n' a = n a + B m a + BM a + \B\ 2 l a , (25) 
m' a = m a + Bl a . (26) 



Let us note that this transformation leaves the electromagnetic tensor F a p given in ( pl| ) invariant; hence, this is the 
gauge subgroup of the little group of l a . 

The rotation part of the subgroup under discussion is simply given by n' a = n a and m' a = Cm a , where C = 
exp(— i<d). This corresponds to a simple rotation by a constant angle in the (e^s, e %)) plane, i.e., 



eft = cos e Q (1) + sin e Q (2) , (27) 
eg } = - sin e a {1) + cos e Q (2) . (28) 

Finally, we note that a mt/Z rotation is an element of the four-parameter group of transformations given by 

l' a = Al a , (29) 
n' a = A^n" + B m a + Bm a + \B\ 2 Al a , (30) 
m ' a = C (m a +B Al a ) . (31) 

A null rotation is the most general transformation of the local null tetrad frame that leaves the spatial direction of 
the null vector l a invariant. 

V. EQUATIONS OF LIGHT GEODESICS AND THEIR SOLUTIONS 

We consider the motion of a light particle ( "photon" ) in the background gravitational field described by the metric 
(^). No back action of the photon on the gravitational field is assumed. Hence, we are allowed to use equations of 
light geodesies ( |l9| ) with l a = dx a /d\ directly applying the metric tensor in question. Let the motion of the photon 
be defined by fixing the mixed initial-boundary conditions 

x(*o)=xo, ^r^ = k ' ( 32 ) 

where |k| 2 = 1 and, henceforth, the spatial components of vectors are denoted by bold letters. These conditions 
define the coordinates x of the photon at the moment of emission of light, i , and its velocity at the infinite past 
and infinite distance from the origin of the spatial coordinates (that is, at the so-called past null infinity denoted by 

j- M)- 

In the underlying inertial frame of the background flat space-time the unperturbed trajectory of the light ray is a 
straight line 

x\t) = x l N (t) = x l + k l (t - t ) , 



where to, x , and (k l ) = k have been defined in equation (|32|). It is convenient to introduce a new independent 
parameter r along the photon's trajectory according to the rule H , 

t = k • xat(<) = t — t + k • x . (33) 



The time to of the light signal's emission corresponds to r = To, where To = k • xo, and t — corresponds to the 
coordinate time t — t* , where 

t* = t - k • x . (34) 

This is the time of the closest approach of the unperturbed trajectory of the photon to the origin of an asymptotically 
fiat harmonic coordinate system. We emphasize that the numerical value of the moment t* is constant for a chosen 
trajectory of light ray and depends only on the space-time coordinates of the point of emission of the photon and the 
point of its observation. Thus, we find the relationship 

T = t-t* , (35) 

which reveals that the differential identity dt — dr is valid and, for this reason, the integration along the light ray's 
path with respect to time t can be always replaced by the integration with respect to t with the corresponding shift 
in the limits of integration. 

Making use of the parameter r, the equation of the unperturbed trajectory of the light ray can be represented as 

x\t) =x i N {r) = fcV + f . (36) 

The constant vector (£ l ) = ^ = kx (xo x k) is called the impact parameter of the unperturbed trajectory of the light 
ray with respect to the origin of the coordinates, d = |£| is the length of the impact parameter. We note that the 
vector £ is transverse to the vector k and directed from the origin of the coordinate system towards the point of the 
closest approach of the unperturbed path of the light ray to the origin as depicted in Figure [l]. 

The equations of light geodesies can be expressed in the first post-Minkowskian approximation as follows (for more 
details see paper (j)) 



£<(t) = ^k a k dih al3 (r^)-d T 



(37) 



where an overdot denotes differentiation with respect to time, d T = d/dr, di = P^d/dtf , k a — (1, k l ), k a = (—1, fej), 
k % = hi, Pij = Sij — kikj is the operator of projection onto the plane orthogonal to the vector k, and h a,3 (T,£) is 
simply h a/3 (t,x) with t = t + t* and x = kr + £. For a given null ray, all quantities on the right side of equation 
( |37| ) depend on the running parameter t and the parameter £ which is assumed to be constant. Hence, equation (|37j) 
should be considered as an ordinary second-order differential equation in the variable r. 

Perturbations of the trajectory of the photon are found by straightforward integration of the equations of light 
geodesies (^). Performing the calculations we find 

x\t) =k i +E i (r) , (38) 
x i (T)=x i N (r)+E i (r)-E i (r ), (39) 
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where r = t — t* and tq — to — t* correspond, respectively, to the moments of observation and emission of the photon. 
The functions 2 l (r) and S*(r) are given as follows 

B*(t) = hak^B^ir) - k a h m (r) - -fc^r) + ^k%k q hP«{T) , (40) 

5*(r) = ^kfAD^ir) - k a B m (r) - ^B oq (t) + ^k%k q B^{r) , (41) 
where it is implicitly assumed that h a ^(— oo) = 0, and the integrals B Q/3 (r) and D a ^(r) are given by 

B Q "W = <W + < (r), <V) = f htffoQdv, B a /(r)= f hf(a,£)da, (42) 

J — oo J — OO 

^(r)=<(r) + Df(r), 25#(t) = f B$(a,®da , D?(t) = f B a /{a,i)da. (43) 

J — oo J — oo 

Integrals (^|) and ( ^3|) and their derivatives are calculated in Appendix C by extending a method developed in (t| 
and frTq| . It is important to emphasize, however, that in the case the body moves along a straight line with constant 
velocity and spin, one can obtain the same results by direct computation without employing the methods used in 
Appendix C. 

Equation (|3^) can be used for the formulation of the boundary-value problem for the equation of light geodesies 
where the initial position, x = x(i ), and final position, x = x(i), of the photon are prescribed for finding the solution 
of the light trajectory. This is in contrast to our original boundary-value problem, where the initial position x of the 
photon and the direction of light propagation k given at the past null infinity were specified. All that we need for the 
solution of the new boundary- value problem is the relationship between the unit vector k and the unit vector 

K = -f^, (44) 
|x - x | 

which defines a geometric (coordinate) direction of the light propagation from the observer to the source of light as 
if the space-time were flat. The formulas ( |39|) and (^|) yield 

fc J = -ir-/? l (T,|) + /f(T ,£) , (45) 

where the relativistic corrections to the vector K l are given by 

^ jwnH : ¥° j H, (46) 

I x (t-)- x o| 

We emphasize that the vectors (f3 % ) = /3(r, £) and (/3g) = /3(to, £) are orthogonal to k and are evaluated at the points 
of observation and emission of the photon, respectively. The relationships obtained in this section are used for the 
discussion of observable relativistic effects in the following sections. 
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VI. GRAVITOMAGNETIC EFFECTS IN PULSAR TIMING, ASTROMETRY AND DOPPLER 

TRACKING 

A. Shapiro Time Delay in Binary Pulsars 

We shall give in this paragraph the relativistic time delay formula for the case of the propagation of light through 
the non- stationary gravitational field of arbitrary-moving and rotating body. The total time of propagation of an 
electromagnetic signal from the point Xo to the point x is derived from equations (|3^) and (fll|). First, we use equation 
( |39| ) to express the difference x — xo via the other terms of this equation. Then, we find the total coordinate time of 
propagation of light, t — to from 

t-h= |x-x | + A M (Mo) + AsMo) , (47) 

where |x — Xo| is the usual Euclidean distance between the points of emission, Xo, and observation, x, of the photon, 
Ajvf(i, to) is the Shapiro time delay produced by the gravitoelectric field of a point-like massive body, and As(i, to) is 
the Shapiro time delay produced by the gravitomagnetic field of the spinning source. The term Ajvf(i)io) is discussed 



in detail in [|53|. The new term As(t, to) is given by (cf. equation (C16) of Appendix C) 



A s (t,to) = B s (T)-B s (T ) , (48) 
B s (r) = -k a k,B s (r) = 2 ^ _ = ^ - - ■ - - ■ - , (49) 

Bs(ro) = \k a k,Bf{ro) = 2 i^ = ^ ? k ^ f , (50) 

2 D y/\ - vl [To - v • r )(r - k • r ) 

where the times r = t — t* and tq = to — t* are related to the retarded times s and so via equation (]l2]), r = x — z, 
i"o = xo — zo, ro = I ro I , x = x(i), xo = x(io), z = z(s), Zo = z(so), v = v(s), and vo = v(so). It is worthwhile to 
note that in the approximation where one can neglect relativistic terms in the relationship ( ^5|) between vectors k 
and — K, the expression for r is given by 

r = Dk + x -z(s) , (51) 

where D = |x — xo|. In the case of a binary pulsar, the distance D between the pulsar and the solar system is much 
larger than that between the point of the emission of the radio pulse, Xo and the pulsar or its companion z(s). Hence, 

r-k, = - [kx( ^ z)]2 k, (52) 

where terms of higher order in the ratio |xo — are neglected. 

The result ([i"8| ) can be used, for example, to find out the spin-dependent relativistic correction As to the timing 
formula of binary pulsars. The binary pulsar consists of two bodies — the pulsar itself (index "p") and its companion 
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(index "c"); in what follows we let r p = x — z p (s), r c = x z c (s), ro p = xo — z p (so), and ro c = xo — z c (sq), where 
z p and z c are coordinates of the pulsar and its companion, respectively. According to Jig] , the difference between 
the instants of time s and sq, s and to, an d so and to for binary pulsars is of the order of the time required for light 
to cross the system, that is of the order of a few seconds. Thus, we can expand all quantities depending on time in 
the neighborhood of the instant to and make use of the approximations r p — x — z p (to) = p p , r c = x — z c (tp) = p c , 
r 0p = x - z p (t ) = p 0pl r 0c = x - z c (t ) = p 0c , J p = J 0p , and J c = J Gc . The next approximation used in the 
calculations is 

xo = z p (<o) +kl , (53) 

where X is the distance from the pulsar's center of mass to the point of emission of radio pulses fl54| . One can also 
see that p 0c = R + hX, where R — z p (t ) — z c (t ), i.e. the radius- vector of the pulsar with respect to the companion. 
Taking these approximations into account, the following equalities hold 

p c -kp c = -2i|^k, (54) 
Pp-kp p = . (55) 

Moreover, 

»-k. ft - *-^' -*^Li«g. (56, 
2p c 2p c 

The primary contribution to the function As(t,to) is obtained after expansion of the expressions for B(t) and 



B(tq) in powers of v/c and picking up all velocity-independent terms. Accounting for equations ( p4| ) - ( |56| ) and (D3) 
this procedure gives the time delay correction As to the standard timing formula as follows 

A ^ _ 2J p ■ (k x v p ) + 2J C ■ (k x r e ) _ 2J 0p ■ (k x r 0p ) _ 2.70c • (k x r 0c ) 
r p (r p - k • r p ) r c (r c - k • r c ) ro p (ro p - k • r 0p ) ro c (r 0c - k • r 0c ) 

By means of ( |53"| ) and (^) we find that the first and third terms in (^?j) drop out. Neglecting terms of order X/R we 
can see that p Qc — R, which brings As into the form 

1 1 



A s = 2^ 0c • (k x R) 

p c (p c -k-p c ) i?(i?- k-R) 
Making use of (|5^) transforms expression (|58| ) to the simpler form 



(58) 



= 2^ 0c ■ (k x R) 2 J-qc ■ (K x R) 
S i?(i? + k-R) i?(i?- K-R) ' 1 j 

where the unit vector K is defined in (f44|). 

Formula ( |5^ ) coincides exactly with that obtained on the basis of the post-Newtonian expansion of the metric 
tensor and subsequent integration of light-ray propagation in the static gravitational field of the pulsar companion 
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p5| . In principle, the additional time delay caused by the spin of the companion might be used for testing whether 
the companion is a black hole or not [Q. However, as shown in [p5| , the time delay due to the spin is not separable 
from the delay caused by the bending of light rays in the gravitational field of the companion [|57j . For this reason, the 
delay caused by the spin is not a directly measurable quantity and can not be effectively used for testing the presence 
of the black hole companion of the pulsar jHf| . 

B. Deflection of Light in Gravitational Lenses and by the Solar System 

Let us assume that the observer is at rest at an event with space-time coordinates (i,x). The observed direction s 
to the source of light has been derived in [Q and is given by 

s = K + a + /3-/3 + 7, (60) 



where the unit vector K is given in (44), the relativistic corrections j3 and /3 are defined in (J46J), a. describes the 
overall effect of deflection of the light-ray trajectory in the plane of the sky, and 7 is related to the distortion of 
the local coordinate system of the observer with respect to the underlying global coordinate system used for the 
calculation of the propagation of light rays. 

More precisely, the quantity a. can be expressed as 

a* = a l M + a l s , ot M = -diB M (r) + k a P)h a A f(r) , 4 = ~9,B s (t) + k a P)h a s \r) , (61) 

the quantities /3 = (3 M + f3 s and /3 = /3 0M + (3 0S are defined by ( [46| ) , and 

y=7lf + 7L Ym = -lk n P)h n A j(r) , Ys^-\k n P)h^(r). (62) 

In what follows, we neglect all terms depending on the acceleration of the light-ray-deflecting body and the time 
derivative of its spin. Light-ray deflections represented by ot l M , [3 l M , and caused by the mass-monopole part of the 
stress-energy tensor ([j]) have been calculated in |Te| 1 and will not be given here as our primary interest in the present 
paper is the description of the spin-dependent gravitomagnetic effects. Using equation (E9J) for the function E>s(t) 



and formula (D£), we obtain 

'(1 - v 2 )k a rpS a Pp i3 r3 (1 - k • \)k a r S a ^ ' P^ 



1 k • v 

dMr) = -2 - 7 _ f 
V 1 — V 



(r — v ■ r) 3 (r — k ■ r) (r — v • r) 2 (r — k • r) 2 



(63) 



(r — v • r) 2 (r — k • r) (r — v • r) (r — k • r) 
In addition, making use of formula ( pll ) yields 



WflgV)-! ^ ^f" + '^f , (64, 

2 -y/l — v 2 r — v • r (r — k ■ vy — v 2 r — k • r 
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2 a y/l - vl r - v • r (r - k • r ) 2 ^/l - r - k • r 



2 k a Pij Sq ^ 



(65) 



The light deflection vector, i.e. the vector connecting the undeflected image to the deflected image of the source of 
light, in the plane of the sky is defined in Jig] and for the purely spin-induced part can be calculated from the third 
relation in equation (|6l|). In that relation the the first term is given in ( |63| ) and the second one can be calculated 
from equation (jig), which results in 



Pi,k a h" 3 = -2Vl - v 2 



(1 _ k . V ) - P„V kar ? Sa0 



(66) 



(9 — v • r) 3 (9 — v • r) 3 

The most interesting physical application of the formalism given in the present section is the gravitomagnetic 
deflection of light in gravitational lenses and by bodies in the solar system. In both cases the impact parameter d of 
the light ray is considered to be extremely small in comparison with the distances from the body to the observer and 
the source of light. The gravitational lens approximation allows us to use the following relationships (for more details 
see section 7B in flq]) 



kr 



r + kr = Co , 



(67) 



where C = P l Ax^ — z 3 (s)) and Q — P l j(x 3 — z J (so)) are impact parameters of the light ray with respect to the light- 
deflecting body evaluated at the moments of observation and emission of light. Assuming that |£| = d <C min[r, j"o] 
and |Col = do <ti mm[r, ro], where r — |x — z(s)\ and r$ = |xo — z(sq)| are distances from the body to observer and to 
the source of light, respectively, one can derive the following approximations 



? — k • r = — 

2r 



r - k • r = 2r 



(68) 



along with 



r — v • r = r(l — k • v) , 



r - v • r = r (l + k • v ) 



Hence, the relativistic spin-induced deflection of light a l s in ( |6l| ) can be calculated from 

2 k a r pS a ^ Q k a PijS° 



d i 



d 2 



(69) 



(70) 



where we have neglected all residual terms of 0(d/r). After substituting the expressions given in (D2) and (D3) in 
( f70| ) , one obtains the analytic representation of the relativistic deflection of light valid for the body having arbitrary 
high speed v and constant spin 3 



87^ 
<i 4 



J-(kx() + J.((xv) 



1-7 
jv 2 



(vJ)(kxO-v 



c 



472 
If 



(v x v)' 



(71) 
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In the case of slow motion, the Taylor expansion of (|7lj) with respect to the parameter v/c yields 

i 8j-(kxC)C , Hkxjy 

as = * + cP ' 



(72) 



which exactly coincides with the previously known formula for a stationary rotating gravitational lens (see, e.g., the 
leading terms of equation (6.28) in j58[) derived using a different mathematical technique. One can also recast ( fr^) 
in a simpler gradient form 

4 = 4^f , ^ s = 4(k x jy din d/8? , (73) 

where ?/>s is a gravitomagnetic component of the gravitational lens potential (cf. the second term of equation (153) 
in @). 

C. Gravitational Shift of Frequency and Doppler Tracking 

The special relativistic treatment of the Doppler frequency shift in an inertial system of Cartesian coordinates is 
well-known. It is based on two facts: proper time runs differently for identical clocks moving with different speeds and 
electromagnetic waves propagate along straight lines in such an inertial system in flat space-time (see ]5^ ] and ]6C| ] 
for more details). The general relativistic formulation of the frequency shift in curved space-time is more involved. 

Two definitions of the Doppler shift are used pl[ — in terms of energy (A) and in terms of frequency (B) 

(A) - = (B) ^ = §, (74) 

v Q u%l 0a vq dT 

where Vq and v are emitted and observed electromagnetic frequencies of light; here (To, Uq ) and (T, u a ) are respectively 
the proper time and 4-velocity of the source of light and observer, and l$ a and l a are null 4-vectors of the light ray 
at the points of emission and observation, respectively. Despite the apparent difference in the two definitions, they 
are identical as equalities u a — dx a /dT and l a — dtp/dx a hold. In order to connect various physical quantities at the 
points of emission and observation of light one has to integrate the equations of light propagation. 

The integration of null geodesic equation in the case of space-times possessing symmetries has been known for a 
long time and extensively used in astronomical practice (see, e.g., f32f and references therein). However, interesting 
astronomical phenomena in the propagation of light rays in curved space-time are also caused by small time-dependent 
perturbations of the background geometry. Usually, the first post-Newtonian approximation in the relativistic N-body 
problem with fixed or uniformly moving bodies has been applied in order to consider the effects of the N-body system on 
electromagnetic signals pq] . Unfortunately, this approximation works properly if, and only if, the time of propagation 
of light is much shorter than the characteristic Keplerian time of the N-body problem. An adequate treatment of the 
effects in the propagation of light must account for the retardation in the propagation of gravitational field from the 
light-deflecting body to the point of interaction of the field with the electromagnetic signal. 
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A theory of light propagation and Doppler shift that takes account of such retardation effects has been constructed 
in the first post-Minkowskian approximation by Kopeikin and Schafer Jl5| and Kopeikin J63| . For the sake of brevity, 
we do not reproduce the formalism here and restrict ourselves to the consideration of the case of gravitational lens only. 
Some details of this approximation have been given in the previous section. Making use of cither of the definitions 
( |74| ) , we obtain for the gravitational shift of frequency 

— = (-V + — v + — v ) ■(aM + ocs) , (75) 

V 0/gr V U U ' 

where D = |x — xo| is the distance between the source of light and observer, r = |x — z(s)\ is the distance between 
the lens and observer, ro = |xo — z(so)| is the distance between the lens and the source of light, v — dx(t)/dt is the 
velocity of the source of light, Vq = dxo(io) / dto is the velocity of the observer, and olm and as represent vectors of 
the deflection of light by the gravitational lens given in (|6l| ) . Formula (|7^) can be applied to the processing of Doppler 
tracking data from spacecrafts in deep space. In the case of superior conjunction of such a spacecraft with the Sun 
or a planet, the result shown in ([75]) should be doubled since the light passes the light-deflecting body twice — the 
first time on its way from the emitter to the spacecraft and the second time on its way back to the receiver. The 
Doppler tracking formula after subtracting the special relativistic corrections (for details see (]l5) and Q) assumes 
the following form 

'5v\ _{8GM : 16G(k x J) ■ C\ r _. r \ ^ 8G 

gr 

Only terms proportional to the mass M of the deflector were known previously (see, e.g., [JOl and |34j). With 
the mathematical techniques of the present paper general relativistic corrections due to the intrinsic rotation of the 
gravitational lens can now be calculated. For the light ray grazing the limb of the light-ray-deflecting body, the 
gravitomagnetic Doppler shift due to the body's rotation is smaller than the effect produced by the body's mass by 
terms of order w rot L/c, where w ro t is the body's angular frequency and L is its characteristic radius. In the case of 
the Sun, the effect reaches a magnitude of about 0.8 x 10~ 14 , which is measurable in practice taking into account the 
current stability and accuracy of atomic time and frequency standards (~ 10~ 16 , cf. |65| and |66|| ). For Jupiter, the 
corresponding estimate of the gravitomagnetic Doppler shift is 0.7 x 10~ 15 , which is also, in principle, measurable. 

VII. THE SKROTSKII EFFECT FOR ARBITRARY MOVING POLE-DIPOLE MASSIVE BODIES 

A. Relativistic Description of Polarized Radiation 

The polarization properties of electromagnetic radiation are defined in terms of the electric field measured by an 
observer. Let us start with a general electromagnetic radiation field F a p and define the complex field T a p such that 
F a p = Re(^ r ct ^) and E a — Re(£ Q ), where £ a = T a pu^ is the complex electric field. In the rest frame of an observer 
with 4- velocity u a , the intensity and polarization properties of the radiation are describable in terms of the tensor 



SUM LbU Jtx J -c, \ / r r \ SC* / r r \ 

j ,„ [ W + c¥ J ( V - D V D Vo ) - C -^V-D V - D V0 ) - {kxJ) - (?6) 
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Ja/3 — < £a£/3 > ! 



(77) 



where the angular brackets represent an ensemble average and Japu 13 — 0. The electromagnetic Stokes parameters 
are defined with respect to two of the four vectors of the tetrad e a my introduced in (p2|), as follows (cf. |67) and [|68|) 



5*0 = JaP 

51 = J a /3 

52 = Ja/3 

53 = iJ a /3 



(1) (1) + (2) (2) 

P a p@ _ p a pP 
e (l) e (l) e (2) e (2) 

P a pi 3 _|_ P a P P 
C (2) + e (2) e (1) 

a (3 _ a (3 
(l) e (2) e (2) e (l) 



(78) 
(79) 
(80) 
(81) 



Using equation (|77|), the Stokes parameters can be expressed in the standard way |f37jl in a linear polarization basis as 



So = <|£ (1) | 2 + |£ (2 )| 2 >, 



S x = < \£ 



(i)l 



\£, 



(2)1 



> 



S*2 = < £(1)£(2) + £(1)£(2) > > 

5 3 = i < £(i)£(2) —£{i)£{2) > 5 



(82) 
(83) 
(84) 
(85) 



where £, 



(n) 



(n) 



for 71 = 1,2. Under the gauge subgroup of the little group of l a , the Stokes parameters remain 



invariant. However, for a constant rotation of angle O in the (e a m , e °(2)) pl ane i S' = So, S{ = Si cos 20 + 52 sin 20, 
S' 2 = —Si sin 20 + S2 cos 20, and S' 3 = S3. This is what would be expected for a spin-1 field. That is, under a duality 
rotation of = 7r/2, one linear polarization state turns into the other. 
For the null field (j2lj) under consideration in this paper, 



£ a = ~w(<I> m a + * m a ) 



(86) 



where u = —l a u a is the constant frequency of the ray measured by the observer with 4-velocity u a . Changing 
from the circular polarization basis ( |S6|) with amplitudes ujQ/2 and uj^/2 to the linear polarization basis £ a — 
^W e \i) + ^(2) e "2)' we nnc ^ tnat ^(1) = + *)/v / 8 and £( 2 ) = _ &)/V8, The variation of the Stokes 

parameters along the ray are essentially given by equation (|23|) , since the frequency u> is simply a constant parameter 
along the ray given by lj — dt/dX = dr/dX. 

The polarization vector P and the degree of polarization P = |P| can be defined in terms of the Stokes parameters 
(So, S) by P = S/So and P = |S|/So, respectively. Any partially polarized wave may be thought of as an incoherent 
superposition of a completely polarized wave with Stokes parameters (PSq,S) and a completely unpolarized wave 
with Stokes parameters (So — PSo, 0), so that (So, S) = (PSo, S) + (So — PSo, 0). For completely polarized waves, 
P describes the surface of the unit sphere introduced by Poincare. The center of Poincare sphere corresponds to 
unpolarized radiation and the interior to partially polarized radiation. Orthogonally polarized waves represent any 
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two conjugate points on the Poincare sphere; in particular, Pi = ±1 and P3 = ±1 represent orthogonally polarized 
waves corresponding to the linear and circular polarization bases, respectively. 

Any stationary or time-dependent axisymmetric gravitational field in general causes a relativistic effect of the 
rotation of the polarization plane of electromagnetic waves. This effect was first discussed by Skrotskii ( |2j| and 
later by many other researches (see, e.g., |H| an( l @1 an d references therein). We generalize the results of previous 
authors to the case of spinning bodies that can move arbitrarily fast. 

B. Reference Tetrad Field 

The rotation of the polarization plane is not conceivable without an unambigious definition of a local reference 
frame (tetrad) constructed along the null geodesic. The null tetrad frame based on four null vectors (P, n a , m a ,m a ) 
introduced in section 4 is a particular choice. As discussed in section 4, this null tetrad is intimately connected 
with the local frame (e a r \, e °\r\-> e °(2)i e< (3))> wnere the vectors e" ^ and e"^ are defined in ( p2|) and the spacelike 
vectors e"-^ and e a , 2 ) are directly related to the polarization of the electromagnetic wave. Each vector of the tetrad 
(e° pe° 1 j,e" 2 j,e" 3 j) depends upon time and is parallel transported along the null geodesic defined by its tangent 
vector l a . To characterize the variation of e"^ along the ray, it is necessary to have access to a fiducial field of tetrad 
frames for reference purposes. To this end, let us choose the reference frame based on the set of static observers 
in the background space-time. Then, at the past null infinity, where according to our assumption the space-time is 
asympotically flat, one has 

ef 0) (-oo) = (1,0,0,0), ef 1) (-oo) = (0,a 1 ,a 2 ,a 3 ), ef 2) (-oo) - (0, b\ b\ b 3 ) , ef 3) (-oo) = (0, k\ k\ k 3 ) . (87) 

Here the spatial vectors a = (a 1 , a 2 , a 3 ), b = (b 1 , b 2 , 6 3 ), and k = (k 1 , k 2 , k 3 ) are orthonormal in the Euclidean sense 
and the vector k defines the spatial direction of propagation of the light ray (see equation (|32|)). At the same time, 
the four vectors ( |87| ) also form a basis e^j = e^d / dx a of the global harmonic coordinate system at each point of 
space-time, where e"^ = e"^(— 00). However, it is worth noting that this coordinate basis is not orthonormal at 
an arbitrary point in space-time. Nevertheless, it is possible to construct an orthonormal basis at each point 

(ft) ~ is -7 c (0) 



by making use of a linear transformation <jJ°iq\ — A" eTm such that the transformation matrix is given in the linear 



approximation by 



A° = l + ^oo, A°. = ifcoi, A l = -h 0z , .V ,>', . 1 ssi 



A local orthonormal basis is then defined by 



(11 = M + ^00 , , , j , = (hoja 3 , a 1 - ^hiji 



(89) 
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By definition, the tetrad frame e"^ is parallel transported along the ray. The propagation equations for these vectors 
are thus obtained by applying the operator D\ of the parallel transport (see (^)). Hence, 

de a 

-^T + r £Aw= ' ( 9 °) 



where A is an affine parameter along the light ray. Using the definition of the Christoffel symbols (Bl) and changing 
over to the variable r with dx a /dr = k a + 0(h), one can recast equations ( p0[ ) in the form 

Equation ( |9l"l) is the main equation for the discussion of the Skrotskii effect. 

C. Skrotskii Effect 

We have chosen the null tetrad frame along the ray such that at t = — oo the tetrad has the property that 
e 00 ) = 0. It follows that in general e ^ = 0(h). The propagation equation for the spatial components e 1 ^ is 
therefore given by 

Furthermore, we are interested only in solving equations (|9l]) for the vectors e?U and e^x that are used in the 
description of the polarization of light; hence 

4- U M + \ h ^(S) + ^ (n P = . (« = 1,2) , (93) 



where we have defined the quantity Q l as 

- -Eijldj (^hlakA • (94) 

Therefore, eVx can be obtained from the integration of equation (p3|). Moreover, ( a e" n j = implies that 

It is worth noting that as a consequence of definition l a = lo(1,x % ) and equation fl38|), one has 1 = u>(k + E). Moreover, 
from the condition l a l a = 0, it follows that k • S = — (l/2)h a pk a k 13 . 

Let us decompose fi* into components that are parallel and perpendicular to the unit vector k l , i.e., 

= (k • fyk 1 + P)Cl j . (96) 

Then, equation (|9^) can be expressed as 

( e \n) + 2% e %)) + ^ ■ "'-'/" !„)" -i--..-/' 1,,,' v 



' eV, + Ihijet.s) + (k ■ O^e^fc' + e^^PK fl« = . (n = 1, 2) (97) 
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Integrating this equation from — oo to r taking into account initial conditions (87) and equalities £ijia 3 k l — —b % and 
EijiVk 1 — a 1 , we obtain 

e 4 (1) = a 1 - ^hija? + (( k-ndajb 1 - e^a 3 P l q f VL q da , (98) 



= & - -hijb 3 - ( k • n da ) a 1 - £ %il VP l q / VL q da . (99) 



e (2) 

To interpret these results properly, let us note that a rotation in the (wV^, ^(o\) plane by an angle <f> at time r 
leads to 

w '(i) = ( a * _ \hijOp) cos0 + (b l - -hijW)sm(j) , (100) 
lo\ 2 ) = -(a 1 - ^hij a 3 ) sin <f) + (b l - ^ h tj V ) cos tf> , (101) 

so that if the angle <fi — O(h), then we have 

uj\ 1} = a 1 - \-hija? + <j> b* , uo\ 2) = b l - ]-h i:j b j - <j> a* . (102) 



Comparing these vectors with equations ( |98| ) and (|99|), we recognize that as vectors em and e( 2 ) propagate along the 
light ray they are rotating with an angle 

4>(t) = ( k • ft da (103) 

J — oo 

about k in the local (up) ,^(2)) plane. Moreover, etn and en) rotate to O(h) toward the direction of light propagation 
k (see the very last terms on the right-hand sides of equations (p8|)and (p9|)). 

We are mostly interested in finding the rotational angle <f> in the plane perpendicular to k. It is worth noting 
that the Euclidean dot product k • ft can be expressed in terms of partial differentiation with respect to the impact 



parameter only. This can be done by making use of equation (C4) and noting that Eij p k 3 k p = 0, so that 

1 



k O = ~k a k l s m d q h afi , (104) 



where the hat over spatial indices denotes the projection onto the plane orthogonal to the propagation of light ray, 
for instance, A 1 = P l jA 3 ■ Hence, the transport equation for the angle <j> assumes the form 



— — k Sipqdqh a fi , (105) 



1 

dr 2 

which is useful for integration. Formula ( |105| ) constitutes a significant generalization of a result that was first discussed 
by Skrotskii ( |2j| and |£0|) and bears his name. 

For a stationary "pole-dipole" source of the gravitational field at rest at the origin of coordinates, equations (|Tc|), 
(pd|), and ( |l6| ) for the metric perturbations imply that 
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hoo = — , hoi = ^ , hij = dij . (10b) 



It follows from (M) that in this case 



„ „ m(k x r) . 
n = B g - 1 3 ; , (107) 



where B ff is the dipolar gravitomagnetic field associated with the source 

'Jxr\ J 



Bg = V X 



KvJ)r-J , (108) 



where f = r/r and J7" = 3 '/ 'J are unit vectors, and J is the magnitude of the angular momentum of the source. 
Thus, k • = k • B g , and, in this way, we recover Skrotskii's original result dcjy/dr = k • B g , which is the natural 
gravitational analog of the Faraday effect in electrodynamics. The Skrotskii effect has a simple physical interpretation 
in terms of the gravitational Larmor theorem fflj| . We note that in the particular case under consideration here, 
the gravitomagnetic field can be expressed as B g = — V(j7" • r/r 3 ) for r ^ 0. Thus for radiation propagating in the 
spacetime exterior to the source, <fi = 3 ■ r/r 3 + constant. In particular, it follows from this result that the net 
angle of rotation of the plane of polarization from — oo to +oo is zero. This conclusion is confirmed later in equation 



(114) as well. 



In the general case, where bodies generating gravitational field are both moving and rotating, the integration of 



equation (105) is also straightforward and is accomplished with the help of a mathematical technique described in 
Appendix B. The result is 

0(r) = 0o + S<P(t) - S(f>(T ) , S<P(t) = ^k a k l e lP gd q B aP (T, £) , (109) 

where </>o is a constant angle defining the orientation of the polarization plane of the electromagnetic wave under 
discussion at the past null infinity, 84>{t) and 8<P(tq) are relativistic rotations of the polarization plane with respect 
to its orientation at infinity at the moments of observation, r, and emission of light, To, respectively, and the function 
B a p is defined in (fi"2j). 

We note that according to the definition ( [i"2"| ) the tensor function B a p consists of two parts the first of which B^f 
relates to the action of the mass monopole field of the particles and the second one Bg describes the action of spin 
dipole fields on the rotation of the polarization plane. Therefore, the angle of rotation can be represented as an 
algebraic sum of two components 

S4> = 64> M + 6</> s , (110) 



where the monopole part can be calculated using equation ( ClCj ) 



hi = \k a k t e^d q B a J = 2m ^-=L| , ^ ( ^ ® r . (Ill) 

2 Vl - v 2 [r - k • r)(r - v • r) 
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The spin part is given by 

Hs = \k a he^d q Bf , (112) 
and has a rather complicated form if written explicitly by making use of partial derivative of function B'g 3 given in 



(p9[). In the gravitational lens approximation when light goes from — oo to +00, equations (111) and (|112j) simplify 
and assume the form 

AGm 

%(+oo) = -^k.(^v), (113) 



. , 4G 
;(+00) = 



(k . (7)+ Ckx€)-(€x^) 



d 2 



. (114) 



Equations (113) and (114) make it evident that in the case of gravitational lensing the integrated effect discussed 
by Skrotskii |2fj] is only due to the translational motion of the lens. There is no contribution to the effect caused by 
spin of the body and proportional to 1/d 2 , where d is the impact parameter of light ray with respect to the body. 
This rather remarkable fact was first noted by Kobzarev & Selivanov |7^] who criticized final conclusions of the paper 
by Skrotskii (25) and others (see (7j| and j74|). We would like to stress, however, that all of the results of the paper 
by Skrotskii are correct, since the final formulas of his paper relate to the situation where light is emitted from 
or near the surface of the rotating body. Such a case obviously is not reduced to that of gravitational lens. Thus, the 
angle of rotation of the polarization plane does not equal to zero as correctly shown by Skrotskii. 

The absence of the Skrotskii effect for radiation propagating from —00 to +00 in the case of a stationary rotating 
body reveals that the coupling of the polarization vector of the photon with the spin of the rotating body can not be 
amplified by the presence of a gravitational lens. 

VIII. THE SKROTSKII EFFECT BY QUADRUPOLAR GRAVITATIONAL WAVES FROM LOCALIZED 

SOURCES 

A. Quadrupolar Gravitational Wave Formalism 

We shall consider in this section the Skrotskii effect associated with the emission of gravitational waves from 
a localized astronomical system like a binary star, supernova explosion, etc. For simplicity we shall restrict our 
considerations to the quadrupole approximation only. The direct way to tackle the problem would be to use the 
Taylor expansion of the formula for the Skrotskii effect given in the previous section. However, it is instructive to 
make use of an approach outlined in |?j that is based on the multipole expansion of radiative gravitational field of a 
localized astronomical source [fr5| . 

To this end, let us consider the propagation of light ray taking place always outside the source with its center of 
mass at rest that emits gravitational waves. The metric perturbation h^ v can be split into a canonical part h ®™', 
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which contains symmetric trace- free (STF) tensors only (for more details on STF tensors see and a gauge part, 

i.e. h uv = ft.™"' + d n w v + d v w n . In the case of mass-monopole, spin-dipole, and mass-quadrupole source moments 
and in the harmonic gauge the canonical part of the metric perturbation is given by 



7 can. 
"■00 



2M 



0„ 



^Lpq (t 



(115) 



■Lean. 
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2d, 



l l3 {t 



(116) 



ican. c ican. 



2 ■• 

~2y (t - r) 



(117) 



Here the mass M and spin J % of the source are constants, whereas its quadrupole moment 2Jy is a function of the 
retarded time t — r. Dependence of the mass and spin on time would be caused by the process of emission of energy 
and angular momentum that are then carried away from the source by gravitational waves. If necessary, the time 
dependence of mass and spin can be treated in the framework of the same calculational scheme as applied in the 



present paper. Moreover, it is assumed in (115) - (117) that the center of mass of the source of gravitational waves 
is at the origin of the coordinate system and does not move. Hence, the radial coordinate r is the distance from the 
center of mass to the field point in space. 

The explicit expressions for the gauge functions u> M relating h^J 1 ' with h uu are chosen in the form [|l5| 



C-%iw(t-r) 



(118) 



1. 



w 



-2V fe 


'lki(t - r)~ 




r 



where we have introduced the definitions of time integrals of the quadrupole moment 

= /" dvlij(v), ^I l3 (t)= [ dv^X i3 {v) 



(119) 



(120) 



Gauge functions (118) and (119) make space-time coordinates satisfy both harmonic and ADM gauge conditions 
simultaneously . The harmonic- ADM coordinates are especially useful for integrating equations of light propagation 
and the description of the motions of free falling source of light and observer. It turns out that the source of light and 
the observer do not experience the influence of gravitational waves in the ADM coordinates and move only under the 
action of the stationary part of the gravitational field created by the mass and spin of the source of the gravitational 
waves (for more details concerning the construction of the harmonic- ADM coordinates see ) . 

The rotation of the polarization plane is described by equation ( |105| ), where the rotation frequency fl is defined in 
(|§). Using expressions (|115[) - (|117|) together with the gauge freedom, the angular velocity is given by 



n = d T 



± jp (t - r) 




k 3 Ij p (t — r) 




r 


r 





- k ■ (V x w) 



(121) 
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where w = (w l ) is given by (|119|), r — y/r 2 + d 2 , and d = |£| is the impact parameter of the light ray. 



Integration of equation (10S) with respect to time results in 

5<t>(r) : 

where B = (B l ) with 



7"5 k,(VxB) ' 



fl4 = Wl l3 {t-r) + 2Vl l3 {t-r) _ 
ry r 



1ki{t - r) 



(122) 



(123) 



and y = T — r. Making use of formula (119) for vf and taking partial derivatives brings equation (122) to the following 
more explicit form 



(r) 



JiX % (kx 



ry 



T ■ T ■ T 
ry r 2 r 



(k x £) V 



(k x ^) l fc^ 



(124) 



where the quadrupolc moments are all evaluated at the retarded time s = t — r. It is straightforward to obtain 
two limiting cases of this formula related to the cases of a gravitational lens and a plane gravitational wave. 



B. Gravitational Lens Approximation 



In the event of gravitational lensing one has 



y 



d 2 



<P_ 

2r 



(125) 



so that 

d 2 d 2 
H/ = -y, t-r = t* - — + ..., (126) 

where t* is the instant of the closest approach of light ray to the center of mass of the gravitational lens. After making 



use of the approximations shown above, equation (124) simplifies and assumes the form 



5(f>(t* 



4± u (r)(kxg)^ 
d 4 



(127) 



One can see that the first nonvanishing contribution to the Skrotskii effect comes from the time derivative of a 
mass quadrupole moment and does not depend on the spin of the source of gravitational waves. It agrees with the 



conclusions of section 7. Formula (127) also shows that if the source of gravitational waves is periodic like a binary 
star system the polarization plane of the electromagnetic wave will experience periodic changes of its orientation with 
a characteristic frequency that is twice that of the source (72| . 



Equation (127) can be derived from the equation (113) as well. Indeed, let us assume that the gravitational lens is 



comprised of N point particles forming a self-gravitating body. For N point particles, equation (113) implies 

N 



HI 



(+oo) = -4]T 



Gm a 



rk • [(£ - x v a ] 



(128) 
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where the point particles are enumerated by index 'a' running from 1 to N, m a is mass of the ath particle, x a and 
v a are coordinates and velocity of the ath particle, = P^x^, and is the impact parameter of the light ray with 
respect to the origin of coordinate system that we assume coincides with the center of mass of the lens, X % , defined 
by the equation 

N 



(129) 



As a consequence of ( |129] ) , we conclude that all of the time derivatives of X 1 vanish identically. We also define the 
spin of the lens and its tensor of inertia as follows 

N N 



J 1 = ^m a (x a x v a ) 



/«(**) = £m a 4(i*)4(i* 



(130) 



a— 1 a— 1 

The spin of the lens is conserved and does not depend on time, while the tensor of inertia is, in general, a function of 
time. 



Let us expand the right-hand side of equation (128) in a power series with respect to the parameter |£ a |/|£| which 
is presumed to be small. We find that 



(+oo) = - 



AG 



N 



N 



k ' £ X ^ m aVa - k • ^ m °- (£<l X V a) 



N 



^^£"U(€«-fl + 0(^) , (131) 



\ a— 1 / a— 1 

where d — \£\. Taking into account definitions of the center of mass, spin, tensor of inertia of the lens, as well as the 
identity 



2[(k x €) • Vo ](£ a • = (k x & tjC -(£xJ) 



(132) 



together with formula (113) we arrive at the result shown in ( |127|) . Comparing expression (127) with the corresponding 
result derived by Kobzarev & Selivanov ( j72|, formulas (11) and (12)), we conclude that the result given by these 
authors |Q is misleading. 



C. Plane Gravitational Wave Approximation 



In the limit of a plane gravitational wave, we assume that the distance D between the observer and source of light 
is much smaller than r and rp, their respective distances from the deflector. The following exact equalities hold: 



d = r sin $ 



y = t — r = r(cos i? — 1) 



(133) 



where z9 is the angle between the directions from the observer to the source of light and the source of gravitational 
waves. Furthermore, we note that the vector corresponding to impact parameter d, can be represented as 



C =r(N i -k i cos'&) , 



(134) 



2G 



where N l = x % jr and |N| = 1. Making asymptotic expansions up to leading terms of order 1/r and 1/ro, and 
neglecting all residual terms of order 1 jr 2 and 1/r 2 ,, would lead to the following result 

= (kxN)^^-^( t -r) = (kxN)^(t-r) 
rv ' r(costf- 1) r(l-costf) v ; 

where the definition of "transverse-traceless" tensor |47j with respect to the direction N is taken into account 

= lij + i (S i:i + NtNj) N p N q l pq - {SipNjNq + S jp NiN q ) l pq , (136) 

and the projection is onto the plane orthogonal to the unit vector N. In terms of the transverse-traceless metric 
perturbation hj^ = 2Tjj T (t — r)/r, the angle of rotation of the polarization plane of electromagnetic wave emitted at 
past null infinity is given by 

In the case of electromagnetic wave emitted at instant to and at the distance r$ from the source of gravitational 
waves and received at instant t at the distance r, the overall angle of rotation of the polarization plane is defined as 
a difference 



(r) tyfo) = ~ 



(k x Nyy hTT (k x N )^ tt ■ 



(138) 



1 — cos $ lJ 1 - cos i? 

where Nq = Xq/tq and is the angle between the vectors k and Nq. 

It is worthwhile to compare the results obtained in the present section to the standard approach used for calculating 
effects caused by plane gravitational waves. Ordinarily in such an approach a plane gravitational wave is assumed to 
be monochromatic with the property that at infinity the space-time remains asymptotically flat. Actually it means 
that one deals with a localized packet of the waves with boundaries asymptotically approaching to plus and minus 
null infinity uniformly. The metric tensor of such a plane gravitational wave in the transverse-traceless (TT) gauge 
can be expressed as Q 

/ioo = h i = = Re{ajj exp(ip a x a )} , (139) 

where the spatial tensor dij is symmetric and traceless {a^h — 0), Pa = ^gr{— l,p) is the propagation 4- vector of the 
gravitational wave with constant frequency u> gr and constant unit spatial propagation vector p, and aijfP — J77[ . 
One can easily see that such a localized plane gravitational wave does not interact at all with an electromagnetic ray 
propagating from minus to plus null infinity. 

Indeed, the integration of the light ray equation (|37| ) leads to integrals involving the exponential function entering 



(139) along the unperturbed light ray path of the form 

lim 



im / e lp « k T dr = 2w6(p a k a ) = 2 7 rw or 1 ( 5(k • p - 1) , (140) 
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where 5 (x) is the Dirac delta function. Now, if the direction of propagation of electromagnetic ray does not coincide 



with that of the gravitational wave (k ^ p) the argument of the delta function in (140) is not zero, i.e. k • p < 1, 
and the result of integration along the light ray trajectory vanishes. On the other hand, if k = p there is again no 
effect since the right side of equation (^) vanishes due to the transversality of hf? revealing that in the case under 
consideration hj^y = hJTp* = 0. This result has been proved by Damour and Esposito-Farese j78| who studied 
the deflection of light and the integrated time delay caused by the time-dependent gravitational field generated by a 
localized material source lying close to the line of sight j7{| - [^0| . 

The absence of interaction of plane monochromatic gravitational waves with electromagnetic rays propagating 
from minus to plus null infinity makes it evident that all relativistic effects taking place in the gravitational lens 
approximation have little to do with gravitational waves — i.e., only the near- zone gravitational field of the lens 
contributes to the overall effects. The other conclusion is that the plane gravitational wave disturbs the propagation 
of electromagnetic signals if and only if the signals go to finite distances - [§2|. However, another limitation has 



to be kept in mind, namely, as it follows for instance from (138) the standard plane gravitational wave approximation 



(139) is applicable only for gravitational waves with dominant wavelengths X gr <C min(r, ro). This remark is especially 
important for having a self-consistent analysis of such intricate problems as the detection of the solar g-modes by 
interferometric gravitational wave detectors p3|, the theoretical prediction of low-frequency pulsar timing noise (see, 
e.g., - [^7), and references therein), and the anisotropy of cosmic microwave background radiation caused by 
primordial gravitational waves. 

APPENDIX A: GEOMETRIC OPTICS APPROXIMATION FROM THE MAXWELL EQUATIONS 

Here we give a brief derivation of the equations of geometric optics approximation, discussed in section 4, from the 
Maxwell equations for electromagnetic waves propagating in empty space-time. The source-free Maxwell equations 
are given by |67|] 

V a F^ + VpF ia + V 7 F Q/3 = , (Al) 
V p F a>3 = , (A2) 



where V Q denotes covariant differentiation |89| . Taking a covariant derivative from equation (Al) and using equation 



(A2), we obtain the covariant wave equation for the electromagnetic field tensor 

□ g F Q/3 + R al3lS F~< s - R ai F J p + R Pl R\ = , (A3) 

where fj s = V a V Q , R a /3~fS is the Riemann curvature tensor, and R a p = R 1 ai p is the Ricci tensor. 

Let us now assume that the electromagnetic tensor F a p has a form shown in equation (|l7|). We introduce a 
dimensionless perturbation parameter e and assume an expansion of the field of the form 
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F, 



a/3 



Re 



(a a /3 + sba/3 + e 2 c af} + ...) exp 



(A4) 



see p8| for a critical examination of this procedure and its underlying physical assumptions. Substituting the expansion 



(A4) into equation (Al), taking into account the definition l a — dtp/dx 01 , and arranging terms with similar powers of 
e lead to the chain of equations 



l a ap 1 + lp a ia + Z 7 a a p — , 
i (V Q apj + V/3 a 7Q + V 7 a a p) = l a + l@ b ia + l 7 b a p , 



(A5) 
(A6) 



etc., where we have assumed the effects of curvature are negligibly small . Similarly, equation (A2) gives a chain of 
equations 

lf3 a a(} = , (A7) 
V/3 a at3 + ilp b af3 = , (A8) 

and so on. 



Equation (A7) implies that the electromagnetic field tensor is orthogonal in the four-dimensional sense to vector 



l a . Contracting equation ( |A5| ) with l a and accounting for ( A7), we find that l a is null, l a l a — 0. Taking the covariant 
derivative of this equality and using the fact that = since l a — V Q </J, one can show that the vector l a obeys 



the null geodesic equation (|l§|). Finally, equation (A3) can be used to show that 

V 7 (a Q(3 P) + FV 7 a al3 = , 



(A9) 



which immediately leads to equation ( po| ) of the propagation law for the electromagnetic field tensor. In this way, one 
can prove the validity of the equations of the geometric optics approximation displayed in section 4. 

APPENDIX B: THE APPROXIMATE EXPRESSIONS FOR THE CHRISTOFFEL SYMBOLS 



Making use of the general definition of the Christoffel symbols 



1 



d a = d/dx a , 



and applying the expansion of the metric tensor (H) result in the approximate post-Minkowskian expressions 



r oo = -7}d t h 0Q (t,x) , 



di = — 

1 dx % 



r° 



y - -- [dih j(t,x.) +djh 0i (t,x) - d t h i3 (t,x.)] , 



r n 



r oj = 9 [9jh i(t,x) - dih j(t,x) + d t hij(t,x)} 



r ip = o I^i^*'*) + 9 p hij(t,x) - dih jp (t,x.)} 



(Bl) 

(B2) 
(B3) 
(B4) 
(B5) 
(B6) 
(B7) 
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These expressions are used for the calculation of the right-hand side of the null geodesic equation (TL 



APPENDIX C: CALCULATION OF INTEGRALS ALONG THE LIGHT-RAY TRAJECTORY 



In order to calculate the integrals ( ff2| ) and (fl3|), it is useful to change in the integrands the time argument, a, to 
the new one, £, defined by the light-cone equation (|l2|), which after substituting the unperturbed light trajectory (|3 
for x can be expressed as follows 



<y + t* =C + l£ + ka-z(C)| , 



(CI) 



where t*, £, and k are considered as parameters such that the orthogonality condition = holds. The differenti- 



ation of equation (CI) yields the partial derivatives of the retarded time with respect to the parameters 

<K /',,/•' d( or 1 



K 

dt* 



P,, 



r — v • r 



13 d& r — v • r ' dk l 
and the relationship between the time differentials along the world-line of the photon 



r — v • r 



da = dC, 



r — v • r 
-kr 



(C2) 



(C3) 



It is worth noting that in the formula (C2) one can write fyr 3 = (k x (r x k)) 



fc 4 (k • r). Moreover, & = Pyg 



since P l] — 8 lJ — k % W and k^ 1 — 0, i.e. £ l has only two independent components. 

Two principal differential identities applied to any smooth function F(t, x) are used in calculations throughout the 
paper, namely, as proved in 



dF(t,x) 



dF(t,x) 
dt 



x=k(t— t )+x 



and 



dF(t, x) 



dt 



dF(T + t*,kT + tj) 9F(T + r,kr + Q 



dF{a + t*M + £) 



(C4) 



(C5) 



i- Fit* 

t=a+t* ; x=ko-+4 UL 

Relationship ( |C4|) allows one to change the order of operations of partial differentiation and substitution for the 
unperturbed light ray trajectory, while ( |C5| ) shows how to change differentiation from the time t to the parameter t* 
making use of the reparametrization of the light-ray trajectory. Parameters £ l and t* do not depend on time and, for 
this reason, derivatives with respect to them can be taken out of the integrals along the light-ray trajectory. Hence, 
applying ( |5i] ) and flC5| ) to the integral p2] ) and accounting for ( |C3| ) we obtain 

l, Qi(a„,0) 

(C6) 



ot J^oo r — k • r 04* J_ 00 r — k • r r — v ■ r 



where the expression 



-k-r = t*+k-z(C)-C 



(C7) 
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is a function of the retarded argument £ and the (constant) parameter if only. Differentiation with respect to the 



parameters t* and in equation (C6) results in 



(r — v ■ r)(r — k • r) 



-oo (r-k-r) 2 



(C8) 



where r" = (r, r') and r Q = (— r, r 1 ). A useful formula is obtained after differentiating with respect to the impact 
parameter. Specifically, we have 



d,B 



a/3 _ 



(r — v • r)(r — k • r) r — k ■ r (r — v ■ r) 2 



4P«r 3 



(1 - w 2 )r 7 £?( a vfl (i _ k ■ v)r 7 ^M 



(r — v • r)(r — k • r) 2 (r — k • r)(r — v • r) 3 (r — v • r) 2 (r — k • r) 2 
For comparison, in the case of pure monopole particles we have Jl5[ 



d l B^(T) = -4m(l - V 2 ) 1 ' 2 2 — 1 ±— 



(C9) 



(CIO) 



r — v • r r — k • r 

As for the integral D'g 3 , one can see from ([IoJ) and (^lj) that in the calculation of the light-ray trajectory we do not 
need D a @ directly but only its partial derivative with respect to the parameter £\ Thus, using the definition of D 
and the expression for B^ given in (C8), one can prove that 



a/3 

S 



diDf(r) 



r — v ■ r (r — k • r) 



-oo ^ k ' r ) 



d( + 8 



-oo (r-k-r) 3 



The remaining integrals in (C8) and ( Cll ) have the following form 

1 - k • v(C) 



1(8) = 



(r — k • r) 



(n = 2,3) 



(Cll) 



(C12) 



where F(£) is a smooth function of the 4-velocity and/or spin of the bodies at the retarded time. The integral (C12) 



can be calculated by making use of a new variable (cf. equation (125)) 



y = k-r-r = C-i*-k-z(C) 



dy = [1 - k • v(C)] d( 



(C13) 



so that the above integral (|C12[ ) can be expressed as 



F(y) dy _ (-i)"+! /•»« pM d_ (1 

dy \y r ' 



n-1 



As) = (-!)*/ ^^-"—^i F(y)^(-^)dy. (C14) 



Integration by parts results in 



I(s) = 



F(y(s)) 



HO dC 



(C15) 



n — 1 (r — k • r) n 1 n — 1 _/_oo (r — k ■ r) n 1 
The last integral in ( |C15|) can be neglected under ordinary circumstances for it includes the acceleration and/or a 
time derivative of the spin of the body. Omitting such terms in accordance with our general approximation scheme, 
we finally arrive at 
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Bs ^ (r- vr)(r-k-r) 1 - k • v (r - k • r) ' (C16) 



r — v • r (r — k • r) 2 1 — k • v (r — k • r) 2 (1 — k • v) 2 (r — k • r) (1 — k • v)(r k • r) 



APPENDIX D: AUXILIARY ALGEBRAIC AND DIFFERENTIAL RELATIONSHIPS 

In this appendix we give several algebraic relationships that can be used for the calculation of observable effects. 
Making use of (||) and definition (||) we have 

S m = 7 (v x Jf , (Dl) 
S ij = 7e ijk J k + ^^(v ■ J)e ijk v k , (D2) 

k a r p S al3 = 7 [J ■ (k x r) + J ■ ((r - kr) x v)] + ^—^{w ■ J)(k x r) • v , (D3) 

ir 

where 7 = (1 - v 2 )" 1 / 2 and e ijk = e oijk [||. 

In addition, one has the following formulas of differentiation with respect to the impact parameter (a dot over 
any quantity denotes differentiation with respect to time t) 

P r3 ~ f) 

d l r^^ , ft = Pfi (D4) 

r — v • r a£? 

v j P ik r k 



d>y = Pij + — , (D5) 

r — v r 



r — v • r 



(D6) 



_1 \ (l-k-v)iy 

— k • r / (r — v • r)(r — k • r) 2 



1 "\ (1 - v 2 + a • r)Pyr J ' , /',,<"' 



\r — v-ry (r — v ■ r)' 3 (r — v • r)^ 

where a = v is the light-ray deflecting body's acceleration and it has been assumed that r = x z(s) = kr + £ — z(s). 
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It is well-known that equation (Al) can be written in terms of usual (non-covariant) derivatives d a = d/dx a . This is due 



to the complete antisymmetry of the equation with respect to the cyclic permutation of all indices. 
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FIG. 1. Astronomical coordinate system used for the calculations. The origin of the coordinate system is at an arbitrary 
point in space. The unperturbed trajectory of a light ray is defined by the unit vector k directed from the source of light 
towards the observer. The impact parameter of the light ray is defined by the vector £ which is orthogonal to k. Gravitating 
bodies with spin move along arbitrary world-lines. 
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